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Abstract 

[~^ I We extend the Polyakov-loop improved Nainbu-Jona-Lasinio (PNJL) model to 2+1 

flavor case to study the chiral and deconfinement transitions of strongly interacting 
•'-H . matter at finite temperature and nonzero chemical potential. The Polyakov-loop, the 

^ • chiral susceptibility of light quarks {u and d) and the strange quark number susceptibil- 



ity as functions of temperature at zero chemical potential are determined and compared 
with the recent results of Lattice QCD simulations. We find that there is always an 
inflection point in the curve of strange quark number susceptibility accompanying the 
appearance of the deconfinement phase, which is consistent with the result of Lattice 
QCD simulations. Predictions for the case at nonzero chemical potential and finite 
temperature are made as well. We give the phase diagram in terms of the chemical 
potential and temperature and find that the critical endpoint (CEP) moves down to 
low temperature and finally disappears with the decrease of the strength of the 't Hooft 
flavor-mixing interaction. 
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1 Introduction 

QCD phase diagram and thermodynamics has been a subject of intense investigation in 
recent years. Lattice QCD simulation is a principal approach to explore the qualitative 
features of strongly interacting matter and make quantitative prediction of its properties. 
Over the years, this formulation has given us a wealth of information about the phase 
diagram and thermodynamics at finite temperature and limited chemical potential (see for 
example Refs. [1-16]. In response to the Lattice QCD simulations, many phenomenological 
approaches in terms of effective degrees of freedom have been developed to give interpretation 
of the available Lattice data and further to make prediction in the region of phase diagram 
that can't be reached by the Lattice QCD. 

A promising ansatz of this sort approach is the Polyakov-loop improved Nambu-Jona- 
Lasinio model (PNJL) [17-20], which combines the two principal non-perturbative features 
of low-energy QCD: confinement and spontaneous chiral symmetry breaking in a simple for- 
malism. The validity of PNJL model has been tested in a series works by confronting the 
PNJL results with the Lattice QCD data [20-24]. It has been reported that the 2 fiavor {u 
and d quarks) PNJL model can reproduce the result that the crossovers for deconfinement 
phase transition and the chiral phase transition almost coincide [19,20]. For finite chem- 
ical potential, further investigations suggest that the thermodynamics and susceptibilities 
obtained in PNJL model are perfectly in agreement with the corresponding Lattice QCD 
data [20-22] . Recently, the impact of Polyakov-loop dynamics on the color superconductiv- 
ity phase transition and on the pion superfiuidity phase transition within PNJL have been 
explored in Refs. [23] and [24], respectively. 

Although the entanglement of the chiral and the Polyakov-loop dynamics turns out to 
be indispensable to understand the nature of QCD phase transitions and thermodynamical 
behavior, the investigations based on this idea up to now are all performed for the two 
fiavor case with small current quark mass so far. Whether the synthesis of Polyakov-loop 
dynamics with the NJL model when including strange quark works well is still unknown. 
Neverthelss, the study of the critical temperature of the 2-1-1 fiavor QCD phase transition by 
the Lattice QCD simulations with physical masses in the continuum limit have recently been 
reported [15,25]. Therefore it is interesting to extend the 2 fiavor PNJL model to the 2-1-1 
fiavor (i.e., including not only u, d quarks, but also s quark) case to compare with the recent 
results of the Lattice simulations. In addition, it is worth investigating the 2-1-1 fiavor QCD 



phase diagram at finite baryon chemical potential in the NJL model with Polyakov-loop. 

One should note that, in the absence of dynamical quarks, the Polyakov-loop expectation 
value can be taken as an order parameter to identify the color confinement and deconfinement 
in pure gauge theories (see for example Refs. [17-20]). With the introduction of dynamical 
quarks with physical masses, the exact Z{Nc) center symmetry for the gauge fields is lost 
and the Polyakov-loop is no longer an exact order parameter for the transition from the 
low temperature, confined phase, to the high temperature, deconfined phase. However, with 
the distribution function of the quark states, the Polyakov-loop expectation value can still 
describe the evolution from color-singlet tri-quark state to color-triplet single quark state. 
Moreover, the results of Lattice QCD [15] show that the Polyakov-loop is also an useful 
quantity to locate the crossover for deconfinement. Then, the terminology "confinement" 
we take here is the statistical suppression of the color-triplet quark propagation, but not 
the dynamical quark confinement. Similarly, the chiral symmetry is explicitly broken by 
the nonvanishing current masses of quarks. However the chiral condensate is also useful to 
indicate the chiral crossover. 

The purpose of this paper is to extend the 2 fiavor PNJL model to 2+1 fiavor case to 
explore the QCD thermodynamics and phase diagram at finite temperature and nonzero 
chemical potential. The focuses are concentrated on the validity of 2+1 fiavor PNJL at 
finite temperature in comparison with Lattice QCD data and its predictions for the case at 
nonzero chemical potential and finite temperature. The paper is organized as follows. In 
Section 2, we extend the PNJL model to the case of 2+1 flavors. In Section 3, we discuss 
the phase transition of strongly interacting matter at finite temperature and zero chemical 
potential by analyzing the Polyakov-loop, the chiral susceptibility of light quarks, the strange 
quark number susceptibility and other related quantities as functions of temperature at zero 
chemical potential and comparing the results with those of Lattice QCD simulations. In 
Section 4, we make predictions for the phase transitions in the case at nonzero chemical 
potential and finite temperature and give a phase diagram of the strongly interacting matter 
in terms of the temperature and chemical potential. Finally, in Section 5, we give a summary 
and conclusion. 



2 The 2+1 Flavor PNJL Model 



Following Ref. [20], we extend the 2 flavor Polyakov-loop improved NJL model to include 
the strange quark and the 2+1 flavor NJL model [26] with a Polyakov-loop can then be given 



as 



8 

X 

a=0 



-K [det/ (^ (1 + 75) i^) + det/ (^ (1 - 75) ^)] - W ($, $ 
where ■?/' = {ipu,'^d,'^s)'^ is the three-flavor quark fleld, 

D^" = d^ - iA^ with A^" = 6:iA^ , A'^ = QA'i'-^ = -iA. 



• T) 



^<Y 



(2) 



Xa are the Gell-Mann matrices in color space, rho = diag(m„, m(i,ms) is the three-flavor 
current quark mass matrix. Throughout this work, we take m„ = rrid = mi, assuming the 
isospin symmetry is reserved on the Lagrangian level, whereas rris is usually different from 
mi, thus the SU{3)f symmetry is explicitly broken. The quark chemical potential matrix fi is 
chosen to proportional to the unit matrix in our work for simplicity, namely, these three flavor 
quarks have identical chemical potential fi. In the above PNJL Lagrangian, the four-point 
interaction term with an effective coupling strength G is U{3)l ® U{3)ji symmetric, where 
tq = \ \lf and Ta, a = 1, . . . , 8 are the eight Gell-Mann matrices in flavor space. The flavor- 
mixing term with coupling strength i^ is a determinant in flavor space, which corresponds 
to the 't Hooft interaction. It breaks the Ua{^) symmetry and leaves SU{3)l (S> SU{3)r 
unbroken. Z//($,<I>*,T) is the Polyakov-loop effective potential. 

The above description shows apparently that the degrees of freedom for temporal gauge 
fleld expressed in a spatially homogeneous background fleld is explicitly included in the 
PNJL model, compared with the conventional NJL model. The Polyakov-loop dynamics 
represented by this background fleld is controlled by the Polyakov-loop effective potential 
W(<I>, $*,T) and the quarks which couple to the Polyakov-loop. This effective potential can 
be expressed in terms of the Polyakov-loop expectation value (or, in other word, the traced 
Polyakov-loop) $ = (Ti cL)/Nc and its conjugate $* = {TtcL'^)/Nc with the Polyakov-loop 
L being a matrix in color space given explicitly by [20] 

r /"^ ^1 

L{x) = V exp i I dr A4 {x, r) = exp [i/3A4] , (3) 



with /9 = 1/T being the inverse of temperature and A4 = iA^. In the so-called Polyakov 
gauge, the Polyakov-loop matrix can be given as a diagonal representation [19]. The coupling 
between Polyakov-loop and quarks is uniquely determined by the covariant derivative D^ in 
the PNJL Lagrangian in Eq. (1). The trace of the Polyakov-loop, $ and its conjugate, $*, 
can be handled with classical field variables in the PNJL. 

Temperature dependent effective potential Z//($, $*, T) is taken to reproduce the thermo- 
dynamical behavior of the Polyakov-loop for the pure gauge case in accordance with Lattice 
QCD predictions, and it has the Z{3) center symmetry like the pure gauge QCD Lagrangian. 
This Z{3) center symmetry is spontaneously broken when temperature is above some critical 
point (To ^ 270 MeV in pure gauge QCD [20]) and the traced Polyakov-loop develops a finite 
value. In the absence of quarks, $ = $* and the Polyakov-loop serves as an order parameter 
for the deconfinement exactly. 

In previous works, two possible forms for the Polyakov-loop effective potential have been 
well developed. One is a polynomial in $ and $* [20], which is denoted by Upoi{^, <I>*,T) in 
this work and another is an improved effective potential in which the higher order polynomial 
terms in $ and $* are replaced by a logarithm [22,23]. We denote this improved effective 
potential by Uimp{^,^*,T). Both effective potentials are taken in the present work to in- 
vestigate whether our results depend on the details of the Polyakov-loop effective potential. 
These two effective potentials have the following forms 

with 

62(T) = ao + ai (-^ ) +a2(-^ ) +a3(-^), (5) 
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and 

Ui^p{<^^^*,T) ^ _la(T)<I>*<I> + 6(T)ln[l-6<I>*<I> + 4(<l>*-V<I>')-3(<I>*<I>)'] , (6) 

with 

a(T) = ao + ai (^) + a, (^^) , 6(T) = h (^) • (7) 

A precise fit of the parameter a, and bi in these two effective potentials has recently been 
performed to reproduce the Lattice QCD data for pure gauge QCD thermodynamics and 
the behavior of the Polyakov-loop as a function of temperature in Refs. [20] and [23] , respec- 
tively. In these two works, Tq = 270 MeV is chosen to be the critical temperature for the 
deconfinement to take place in the pure gauge QCD. 



After performing the mean field approximation for the Lagrangian of the three-flavor 
PNJL in Eq. (1), we obtain the thermodynamical potential density as 
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J (27r) 

+ Tin [l + 3$e-(^"-^)/^ + 3<|.*e-2(i?„-M)/T ^ ^-3(E^-^^)/T^ 
+ Tin [l + 3<|)*e-('^"+^)/^ + 3$e-2('E^"+M)/T _^ ^-3(Eu+^^)m 
+ Tin [1 + 3$e-(^''-^)/^ + 3$*e-2(^<'-^)/^ + e-3(^<'-^)/^] 
+ Tin [l + 3<|.*e"(-^^+'')/^ + 3<|.e-2(^''+^)/^ + e-3(^d+M)/rj 
+ Tin [1 + 3<l>e-(^=-'^)/^ + 3$*e-2(^»-^)/^ + e-3(^-M)/T] 
+T In [1 + 3<l>*e-(^^+^)/^ + 3$e-2(^»+M)/T ^ g-3{i?.+M)/T] | ^ (g) 

where /i is the chemical potential, T is the temperature, 0, is the chiral condensate of quarks 



with flavor i, and Ei = y'jP^+Mf is its corresponding quasiparticle energy, with constituent 
masses for the quark of flavor i 

Mi = m,- AGcj), + 2K(P, 0^. (9) 

As mentioned above, the breaking of the isospin symmetry is neglected throughout this work. 
We have thus (pu = 4>d = 4>i in the absence of isospin chemical potential. 

Minimizing the thermodynamical potential in Eq. (8) with respective to (pi, 0s, $, and 
$*, we obtain a set of equations of motion 

an dvt dvt dvt 

This set of equations can be solved for the fields as functions of temperature T and chemical 
potential fi. It has been shown in Ref. [20] that, in 2 flavor PNJL model, there is $ = $* 
when the quark chemical potential is vanishing (/i = 0). While for /i 7^ 0, $ and $* have 
different values. 

In the NJL sector of this model, five parameters need to be determined. In our present 
work we adopt the parameter set in Ref. [27], mi = 5.5 MeV, rris = 140.7 MeV, GA^ = 
1.835, KA^ = 12.36 and A = 602.3 MeV, which is fixed by fitting m^ = 135.0 MeV, 
TTiK = 497.7 MeV, m^. = 957.8 MeV and U = 92.4 MeV. 



3 Phase Transition in the Case of /i = and T 7^ 

It has been strongly suggested by current Lattice QCD simulations that the transition from 
low temperature hadronic phase to high temperature quark-gluon-plasma (QGP) phase at 
vanishing quark chemical potential is a continuous, non-singular but rapid crossover [16, 28]. 
It has also been demonstrated that because of the non-singularity of the crossover, different 
observables lead to various values of transition temperature (T^) in the 2+1 flavor QCD 
with physical masses both for the light quarks mi and for the strange quark rris even in the 
continuum and thermodynamical limit [15]. In order to determine the critical temperature 
Tc at vanishing quark chemical potential in the 2 + 1 flavor PNJL model, we consider three 
quantities, which were used to locate the transition point in the Lattice QCD simulations in 
Ref. [15], the traced Polyakov-loop (i.e., the Polyakov-loop expectation value. In the follow- 
ing we just denote it as Polyakov-loop for simplicity), the light quark chiral susceptibility 
and the strange quark number susceptibility. 

The chiral susceptibility of the light quark is defined as 

^' = -ir-2- 11 

In order to obtain a dimensionless quantity and renormalize the divergence of the thermo- 
dynamical potential ^2 in Lattice QCD simulations, this quantity was normalized to the 
following one [15]: 

^ = |tta(r)-x,(0)l. (12) 

In the PNJL model we have 
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A(E,,T,$) 
with A (x, T, $) and / (x, T, $) being defined as 

A{x, T, $) = 1 + 3<I)e-^/^ + 3<l>e-2^/^ + e'^^/^ , (14) 

$e-/^ + 2$e-^-/^ + e-3-/^ 

respectively. Here, we do not distinguish $ and $*, because they have the same value at 
/i = 0. 



The strange quark number susceptibility is defined as [29] 
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Figure 1: Left panel: calculated Polyakov-loop $, scaled chiral condensates for light quarks 0; 
and for strange quarks 0^ as functions of temperature at vanishing quark chemical potential 
in the case of considering the polynomial Polyakov-loop effective potential Upoi, the improved 
Polyakov-loop effective potential Uimp, respectively. Right panel: calculated temperature 
dependence of the derivatives d^/dT, d(f)i/dT and d(^sl&T with the two kinds of Polyakov- 
loop potentials. 

The calculated results of the temperature dependence of the Polyakov-loop, the light 
quark chiral condensate and the strange quark chiral condensate at zero quark chemical 



potential with the usual polynomial Polyakov-loop effective potential lApoi and the improved 
Polyakov-loop effective potential Uimp are shown in the left panel of Fig. 1. As we can see, 
when quarks (including the light and the strange) have physical current masses, the chiral 
phase evolution is not a true phase transition but a continuous crossover. At the same time, 
the coupling of the Polyakov-loop to quarks turns the first-order deconfinement transition for 
the pure-gauge QCD into a crossover as the same as that in the two flavor PNJL case [20]. 
To determine the pseudo-transition temperature, we have also calculated the derivatives of 
the Polyakov-loop and the chiral condensates with respect to temperature. The obtained 
results with the two kinds of Polyakov-loop effective potentials are illustrated in the right 
panel of Fig. 1. In order to confront our results to those in the Lattice QCD simulations with 
physical masses for the 2-|-l flavor QCD [15], we have rescaled the parameter Tq in lApoi from 
270 to 200 MeV and in Uimp from 270 to 215 MeV along the way taken in Ref.[20]. After 
such a rescaling, we obersve that both the Polyakov-loop effective potentials give the same 
pseudo-deconfinement transition temperature 175 MeV (see right panel of Fig. 1). This value 
is consistent with the Lattice QCD result Tc{P) = 176(3) (4) MeV [15], where the numbers 
in the parenthesis indicate the errors. 

After rescaling the Tq in the effective potential, we can determine the pseudo-critical point 
Tc for the chiral phase transition with other quantities. From the appearance of the peaks 
of d(f)i/dT and d(j)s/dT shown in the right panel of Fig. 1, one can infer that the pseudo- 
critical temperature of the chiral phase transition for light quarks Tc{l) are almost the same 
as that of strange quark Tc{s). Such a similarity is quite natural due to the flavor-mixing 
interaction in the Lagrangian shown in Eq. (1). In more details, the critical temperatures 
obtained from the polynomial effective potential are Tc(/) = 212 MeV and Tc{s) = 210 MeV, 
while the corresponding temperatures from the improved effective potential are 226 MeV and 
223 MeV. These values are all relatively larger than the recent Lattice QCD results, which 
are T^{xi) = 192(7)(4) MeV in Ref. [25], 151(3)(3) MeV in Ref. [15] and 169(12)(4) MeV 
in Ref. [29]. This feature is the same as the standard NJL model which also has relatively 
larger chiral transition temperature [30]. 

Comparing the upper panel with the lower panel of Fig. 1, one can find that the Polyakov- 
loop corresponding to the improved effective potential changes more rapidly near the Tc{P) 
than that relevant to the polynomial effective potential. The peak of d^/dT in the lower- 
right panel of Fig. 1 is much narrower and higher, which induces two little bumps in the curves 
of d(f)i/dT and d(f)s/dT at Tc{P), while there is no such bump in the curves in the upper-right 
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panel of Fig. 1 for the polynomial effective potential. These phenomena indicate that, for 
the improved effective potential, the deconfinement transition has an obvious influence on 
the chiral condensate. 

One may have a question about why the peak of the temperature derivative of the 
Polyakov-loop $ obtained in the improved effective potential is much narrower than that 
given by the polynomial effective potential. To explore the underlying physics, we recall that, 
once the dynamical quarks with physical masses are introduced, the Z(3) center symmetry 
is explicitly broken, so that deconfinement phase transition changes from a first order one in 
pure gauge theories to a crossover. It is then natural that the temperature derivative of the 
Polyakov-loop smears to a peak with finite height and finite width. The calculated results 
shown in Fig. 1 manifest that both the polynomial effective potential and the improved ef- 
fective potential give definitely the temperature derivative of the Polyakov-loop peaks with 
finite width and finite height. It indicates that both the effective potentials represent the 
dynamics correctly in some sense. As for that the width given by the improved effective 
potential is narrower, we should note that recent Lattice QCD results show that the nor- 
malized pressure for full QCD with dynamical fermions looks the same as that in the pure 
gauge theories [28], which indicates that it may be the dynamics of gluons that drives the 
QCD phase transition, but not that of the fermions [31]. And the crossover of the phase 
evolution is a quite rapid one [16,28]. The improved effective potential replaces the higher 
order terms in $ and $* in the polynomial effective potential by the logarithm of the Jacobi 
determinant which results from integrating out the six non-diagonal Lie algebra directions 
while keeping the two diagonal ones to represent $ [19,22,23]. It means that, the improved 
effective potential describes the dynamics of the Polyakov-loop (i.e., the gluon dynamics) 
more appropriately and correlates to the effect of the dynamical quarks more slightly. The 
narrower width of the deconfinement crossover with the improved effective potential can thus 
be attributed to that the improved Polyakov-loop effective potential represents better the 
gluon dynamics. 

The calculated results of the temperature dependence of the chiral susceptibility of light 
quarks are shown in Fig. 2. The pseudo-critical point Tc{xi) determined from the peak 
of the chiral susceptibility is 208 MeV for Upoi and 213 MeV for Uimp. These values do 
not exactly coincide with the results of Tc(/)s shown in Fig. 1, which reflects the nature of 
the crossover of the phases. However, the discrepancies of the pseudo-critical points from 
different susceptibilities are within 10 MeV. 
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Figure 2: Calculated results of the renormalized chiral susceptibility mi^^xi/T'^ ^-s a func- 
tion of temperature at zero quark chemical potential for the two Polyakov-loop effective 
potentials. 

Fig. 3 shows the calculated temperature dependence of the strange quark number suscep- 
tibility. First of all, we discuss the case including the flavor-mixing effects with i^A^ = 12.36 
in the Lagrangian in Eq. (1). From the calculated results shown as the thinner curves in 
Fig. 3, one can notice that the Xs/T"^ increases monotonically with temperature. Ref. [15] has 
shown that the pseudo-critical point Tc{xs) can be defined as the inflection point of a suscep- 
tibility curve, viz. the peak of the curve of d{xs/T^)/dT. Our calculated results illustrated 
in the right panel of Fig. 3 manifest that the location of this pseudo-critical point is affected 
by both the deconfinement and the chiral crossovers. Because of the common influences by 
these two crossovers, we find Tc{xs) = 203 MeV for the polynomial potential, which is in the 
between of the deconfinement pseudo-critical point Tc{P) = 175 MeV and the chiral critical 
point Tc(/) = 212 MeV, Tf.{s) = 210 MeV. This phenomenon is more obvious in the result of 
the improved effective potential, where a sharp peak appears at the deconfinement transition 
point Tc{P) = 175 MeV on the curve of d{xs/T'^)/dT. It manifests evidently Tc{xs) = Tc{P). 
In addition, the peak is followed by a little platform coming from the contribution of the 
chiral crossover for light quarks due to the flavor-mixing effects. This phenomenon can be 
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Figure 3: Left panel: calculated results of the scaled strange quark number susceptibility 
Xs/T^ as a function of temperature (thinner curves are for the case with KA^ = 12.36 and 
thicker ones for KA^ = 0). Right panel: plots of the calculated d{xs/T^)/dT. 

understood as the follows. Being not color-singlets, there are almost no free strange quarks 
before the deconfinement crossover takes place. As the temperature increases, the crossover 
for the deconfinement occurs and the strange quark number susceptibility develops to finite 
values abruptly. In general case, the quasiparticle energy of strange quarks becomes smaller 
due to the partial reduction of the condensate of strange quarks which is along with the light 
quarks chiral crossover because of the flavor-mixing interactions. It leads it easier to excite 
quarks, and in turn to increase the value of the quark number susceptibility. However, since 
the improved effective potential represents the gluon dynamics better and the deconfinement 
crossover governed by the gluon dynamics occurs more abruptly as mentioned above, the 
increasing rate of the strange quark number susceptibility with respect to temperature in- 
duced by the deconfinement crossover dominates over that induced by the reduction of the 
quark condensate at the deconfinement pseudo-critical temperature. As a consequence, the 
pseudo-critical temperature given by the strange quark number susceptibility coincides with 
that for the deconfinement crossover. As the temperature increases further, the reduction of 
the quasiparticle energy of the quark begins to play the role, it generates then a platform in 
the curves of d{xs/T'^)/dT. In order to study the influence of the deconfinement crossover 
on the Xs further we turn off the flavor-mixing interaction by setting KA^ = 0. The obtained 
results are illustrated in the thicker curves in Fig. 3. We find again that Xs/T"^ gradually 
develops a finite value with the appearance of the deconfinement phase. Since the Xs/T"^ is 
no longer infiuenced by the light quark chiral crossover any more, the difference between the 
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pseudo-critical temperatures given by the two effective potentials is much smaller than that 
with the ffavor-mixing being taken into account. Furthermore, the curve of d{xs/T'^)/dT for 
the polynomial potential moves down to lower temperature obviously and there are not any 
platforms on the curve of d{xs/T'^)/dT for the improved effective potential. Therefore, we 
could qualitatively understand the results of Lattice QCD simulations in Ref. [15] that the 
pseudo-critical point determined from the strange quark number susceptibility is quite close 
to that determined from the Polyakov loop. Nevertheless, it is not obvious that the quark 
number susceptibility has any relations with the chiral crossover for light quarks, which 
maybe indicates that the flavor-mixing interactions are relatively weak at high temperature. 
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Figure 4: Calculated results of the scaled strange quark chiral condensate as a function of 
temperature with different values of the current mass of strange quark for the improved effec- 
tive potential (left panel) and the temperature dependence of the corresponding Polyakov- 
loop $ at several strange quark masses (right panel). 

In order to investigate the influences of the current mass of strange quark on the proper- 
ties of chiral and deconfinement crossover in the case of including the flavor-mixing effects, 
we depict the scaled strange quark condensate and the Polyakov-loop as functions of tem- 
perature with different values of strange quark current mass in Fig. 4, using the improved 
effective potential Uimp- Obviously, one can find from the figure that the chiral crossover is 
quite sensitive to the variation of the current mass of strange quark. The chiral transition 
temperature Tc{s) defined above decreases from about 225 MeV to 188 MeV for the improved 
effective potential Uimp when the value of rus is reduced from 140.7 MeV to 5.5 MeV. The 
same thing occurs for the polynomial effective potential Upoi- Furthermore, the region of 
crossover for the chiral phase evolution takes place in a more narrow range of temperature. 
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Nevertheless, as for the deconfinement transition, differences between the physical case with 
rUs = 140.7 MeV and the hypothetical case with rris = 5.5 MeV are very small. Conse- 
quently, the difference between the chiral pseudo-critical point Tc{s) and the deconfinement 
pseudo-critical point Tc{P) becomes smaller when the current mass of the strange quark 
decreases. 
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Figure 5: The same as Fig. 2 but with rris = mi = 5.5 MeV. 

Fig. 5 shows the renormalized chiral susceptibility rrii'^Axi/T'^ as a function of tempera- 
ture for the case of m^ = mi = 5.5 MeV, viz. neglecting the flavor difference. As one expects, 
the curves in this figure are much narrower and sharper than those in Fig. 2. The chiral 
pseudo-critical temperature Tc{xi) is 182 MeV for the Upoi and 191 MeV for the Uimp, which 
are 26 MeV and 22 MeV smaller than their corresponding values for the case with physical 
strange quark mass, respectively. 

We have also calculated the temperature dependence of the strange quark number suscep- 
tibility Xs/T'^ with iTis = mi = 5.5 MeV. The obtained results are shown in Fig. 6. From the 
figure we can recognize not surprisingly, Tc{xs) = 179 MeV for the polynomial effective po- 
tential, which is consistent with the chiral pseudo-critical temperature Tc{xi) = 182 MeV and 
the deconfinement critical temperature Tc{P) = 178 MeV. As for the improved effective po- 
tential, there are two peaks located at 169 MeV and 186 MeV in the curve of the temperature 
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Figure 6: The same as Fig. 3 but with rris = mi = 5.5 MeV and i^A^ = 12.36. 
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Figure 7: Calculated derivatives of the strange quark number susceptibility versus the 
temperature in the cases: mi = 5.5 MeV, m,s = 140.7 MeV; m,i = irts = 5.5 MeV; and 
m,i = m,s = 3.0 MeV. 

derivative. We have then the pseudo-critical temperature in the case of m,s = m,i = 5.5 MeV, 
T,(xs) = 169 MeV, 186 MeV, which correspond to the Tc{P) = 168 MeV, Tdxi) = 188 MeV, 
respectively. Fig. 7 shows the temperature derivative of the strange quark number suscep- 
tibility for the improved effective potential with several values of current mass of strange 
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quarks. One can observe from the figure that the peak related to the chiral critical point 
grows drastically as the current mass of strange quarks decreases from 140.7 MeV to 3.0 MeV. 
One should note that in order to confront our results to those in Lattice QCD simulations 
in [15], we have rescaled the parameter Tq in lApoi from 270 to 200 MeV and in lAimp from 
270 to 215 MeV to fix the pseudo-deconfinement transition temperature Tc{P) = 175 MeV. 
However, recent Lattice QCD simulations [25,29] show that there are some uncertainties in 
this pseudo-critical temperature for the 2+1 fiavors with physical quark masses, ranging from 
150 MeV to 190 MeV. We repeat the calculations above with different values of Tq, and find 
that the hierarchy in the pseudo-critical temperatures also exists for the cases with physical 
masses. For example, if we do not rescale the value of Tq and keep it being 270 MeV as in 
the pure gauge field case, we find T(.{P) = 212 MeV and Tc{l) = 242 MeV for the improved 
effective potential. The difference between these two values are obvious. For the polynomial 
effective potential this difference becomes smaller but is also obvious with Tc{P) = 227 MeV 
and Tc{l) = 240 MeV. One interesting thing is that when the value of rris is reduced from 
140.7 MeV to 5.5 MeV and Tq is kept at 270 MeV as well, we find that the deconfinement and 
chiral pseudo-transitions are almost coincident, with Tc{P) = 205 MeV, Tc{l) = 211 MeV for 
the improved Polyakov-loop effective potential Uimp and Tc{P) = 210 MeV, Tc{l) = 211 MeV 
for the polynomial effective potential Upoi- 

4 Phase Transition in the Case of /i 7^ and T ^ 

In the above sections, we have mentioned that in the presence of quark chemical potential fi, 

the (traced) Polyakov-loop $ and its conjugation $* satisfying Eqs. (10) are different from 

each other. Viz. 

cb* — cf) 

A$ = (18) 

2 ^ ' 

develops a finite value. It has also been demonstrated that this difference originates from 
the sign problem of the fermion determinant at finite density, which is unavoidable not 
only in Lattice QCD simulations but also in the mean-field approximation [32]. However, it 
was shown that such a difference between $ and $* at finite density may be not of major 
qualitative importance in determining the phase diagram [20,33]. Therefore, we use the 
average value of $ and $* to indicate the pseudo-deconfinement phase transition at finite 
density, viz. 

_ (T)* _|_ (T) 

$ = ^— . (19) 
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Solving the four equations in Eqs. (10), we can obtain $ and $* as functions of temperature 
T and quark chemical potential fi, and in turn the A$ and $ defined above. In the following 
discussions, we are only concerned about quantities A$ and $, therefore, $ is just denoted 
by $ from now on. 
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Figure 8: Calculated results of the scaled light quark chiral condensate (left panel) and its 
derivative with respect to temperature (right panel) as functions of temperature for several 
values of the quark chemical potential. 



In this section only the improved Polyakov-loop effective potential is implemented and 
To in this potential is also chosen to be 215 MeV. We investigate the cases with physical 
quark masses {mi = 5.5 MeV, rus = 140.7 MeV) and including the flavor-mixing effect 
(i^A^ = 12.36), and keep other model parameters unchanged (except explicit explanations). 
Fig. 8 shows the calculated temperature dependence of the light quark condensate and its 
derivative at several quark chemical potentials. From the figure, one can recognize easily 
that, with the increase of chemical potential, the pseudo-critical temperature of the chiral 
transition determined from light quarks decreases and the crossover becomes more and more 
narrow and abrupt ( as shown in the right panel of Fig. 8) and eventually evolves to a first- 
order transition at the critical endpoint (CEP) {Tcep = 128 MeV, fJ^f^j^p = 308 MeV). From 
the right panel of Fig. 8, one can also notice that the difference between the chiral pseudo- 
critical temperature and that for deconfinement becomes small with the increase of /i, which 
means that the chiral phase boundary drops more rapidly than that of the confinement 
phase does with the increase of chemical potential fi. This phenomenon is consistent with 
the results shown in Figs. 16 and 17 of Ref. [33] for the two-fiavor PNJL model. 

We have also calculated the temperature and chemical potential dependence of the scaled 
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Figure 9: The same as Fig. 8, but for the condensate of strange quarks. 



strange quark condensate and its derivative d(j)s/dT. The obtained resuhs in the case of 
that the chemical potential of strange quarks takes the same values as for the light quark 
condensate are illustrated in Fig. 9. The figure shows evidently that the curves of d(j)s/dT in 
the right panel of Fig. 9 have peaks at the same temperature as the curves of d(f)i/dT in the 
right panel of Fig. 8. This is reasonable because 0^ is influenced by 0/ through the 't Hooft 
flavor-mixing interactions. However, higher temperature is needed to reduce the value of 0^ 
to approach zero since the strange quark has much larger current mass. Therefore, we see 
that there is a broad peak after the restoration of the chiral symmetry for light quarks on 
the curve of d(j)s/dT with /i = 350 MeV in the right panel of Fig. 9, which corresponds to 
the further reduction of the condensate for strange quarks. It should be noticed that there 
are some uncertainties in the location of the broad peak corresponding to fi = 350 MeV 
because of the high chemical potential and temperature, since in these regions, the fermionic 
distribution function is nonvanishing at momentum beyond the cutoff of the model, which 
indicates that those regions are out of the scope of NJL-like models. 

Fig. 10 represents the calculated Polyakov-loop as a function of temperature at finite 
quark chemical potential. The figure manifests evidently that the rapid crossover for the 
deconfinement transition at /i = is smeared by nonzero /i and the crossover gets much 
flatter and milder. This result is qualitatively consistent with the mean-field result shown 
in Fig. 10 of Ref. [32]. Fig. 11 shows the A<l> at finite quark chemical potential. We have 
magnified the value of A$ by ten times in order to compare it to the $ conveniently. As 
one observes that the value of A$ increases with the temperature below the deconfinement 
pseudo-critical temperature Tc{P), while above Tc{P), it decreases and approaches to zero. 
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Figure 10: Calculated results of the Polyakov-loop $ (left panel) and d^/dT (right panel) 
as functions of temperature at several values of quark chemical potentials. 

To discuss the chemical potential dependence of the phase evolution in detail, we show 
the calculated quark chemical potential dependence of chiral condensates 0; and 0^ at T = 
100 MeV in Fig. 12. Since the temperature T = 100 MeV is below the temperature of CEP 
(128 MeV), the chiral phase transition for light quarks is in first-order, which corresponds 
to the divergence of the curve of d(f)i/dfi in the right panel of Fig. 12. Furthermore, besides 
the influence of the chiral restoration for light quarks due to the flavor-mixing, exhibiting a 
discontinuity in the curves of (pi and (ps, there is a smooth crossover for further reduction of 
the strange quark condensate at higher values of the quark chemical potential (the second 
peak in the curve of dcps/d^ is located at about 450 MeV). We also study the influences 
of the 't Hooft flavor-mixing strength K on the strange quark chiral transition with finite 
quark chemical potential and finite temperature. The calculated results at temperature 
T = 100 MeV are illustrated in Fig. 13. The figure shows evidently that, as the flavor- 
mixing strength becomes weaker, the temperature of CEP becomes lower. Consequently, the 
flrst-order chiral phase transition of light quarks displaying a discontinuity in the curve of 
(f)s/4>s{T = 0) in Fig. 13 evolves to a continuous crossover. While at higher chemical potential 
for further reduction of the strange quark chiral condensate, since the chiral symmetry for 
light quarks has been restored, the flavor-mixing interactions that couple different flavors of 
condensates can be neglected. Therefore, these three curves in Fig. 13 coincide with each 
other at large values of the chemical potential. 

We also investigate the quark chemical potential dependence of the Polyakov-loop $ and 
A$. The calculated results at several temperatures are illustrated in Fig. 14. Such an 
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Figure 11: Calculated results of the Polyakov-loop $ (thin curves) and scaled A$ (thick 
curves) as functions of temperature at several values of the chemical potential. 

investigation is extremely interesting if the chemical potential is beyond the scope of the 
current Lattice QCD simulations. It is easy to infer from the figure that the Polyakov-loop 
is obviously suppressed at small T and fi. It indicates a tendency of confinement. This 
tendency is reversed with the increase of T or /i. For example, the value of $ corresponding 
to T = 180 MeV (the thin dotted curve) is much larger than those corresponding to other 
two values of the temperature at low /i, since T = 180 MeV is higher than the deconfinement 
pseudo-critical temperature at zero chemical potential {Tc{P) = 175 MeV). The amplified 
A$ is also shown in Fig. 14, which is much smaller compared with $ in the deconfinement 
phase. 

As a finality, we present the phase diagrams of the 2-1-1 flavor PNJL in Fig. 15 and 
discuss their dependence on the flavor-mixing interaction strengthes (left panel) and the 
current masses of quarks (right panel) for the improved Polyakov-loop effective potential. 
The curves for crossover are determined by the location of the peaks of the derivative of 
"quasi" order parameters with respect to T or /i. As for the left panel of Fig. 15, the 
solid, dotted curves indicate the crossover, the first-order chiral transition for light quarks, 
respectively, separated by the critical endpoint (CEP). The dashed curves correspond to the 
crossover for the deconfinement and we only depict them in the range of /i = ~ 350 MeV. 
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Figure 12: Left panel: calculated quark chemical potential dependence of the scaled chiral 
condensates 0/ and 0^ at temperature T = 100 MeV. Right panel: plots of d(f)i/dfi and 
d(f)s/dfi with repect to the chemical potential. 
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Figure 13: Calculated strange quark condensate as a function of chemical potential with 
different values of KA^ at T = 100 MeV. 

The dash-dotted curves indicate the "further chiral crossover" for strange quarks and we 
only depict them in the temperature range T G [0 , 120] MeV, because NJL-type models are 
problematic in the range of large n and high T as discussed above. One can observe from the 
figure that, as the strength of flavor-mixing interactions becomes weaker, the CEP moves 
down to lower temperature. We also find that, when the flavor-mixing interaction strength K 
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Figure 14: Calculated Polyakov-loop $ (thin curves) and the scaled A<l> (thick curves) as 
functions of quark chemical potential at several values of temperature. 

approaches to zero, the CEP moves toward yU-axis gradually and finally disappears from the 
phase diagram. It should be noted that two dash-dotted curves corresponding to two values 
of KA^ coincides with each other because of the reason mentioned above. The same thing 
occurs for the deconfinement crossover outside the chiral phase boundary of light quarks. 

Furthermore, the dependence of the phase diagram on the current mass of quarks has 
been investigated as well. We neglect differences between the light quarks and the strange 
quark, and study two cases with mi = rUs = 5.5 MeV, mi = m^ = 0, respectively. Their 
corresponding phase diagrams are depicted in the right panel of Fig. 15. Here, i^A^ = 12.36 
is chosen. Since the current mass of strange quark is the same as those of light quarks, the 
strange quark chiral transition coincides with those of light quarks and there is no dash- 
dotted line in the phase diagram, when the current mass of strange quark is reduced from 
140.7 MeV to 5.5 MeV, both the critical value of /i at T = and that of T at ;U = decrease, 
and the CEP moves toward to the T-axis. If the values of both m,i and m,s are further 
reduced, for example in the chiral limit, we find that the CEP disappears from the phase 
diagram and there is only first order chiral transition in the whole range of chemical potential 
and temperature, which is illustrated by the thinner dotted line in the right panel of Fig. 15. 
This result is quite consistent with the Pisarski-Wilczek argument [34] that the order of the 
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Figure 15: Left panel: calculated phase diagram in terms of temperature and quark chemical 
potential of the 2+1 flavor PNJL with physical masses, viz. mi = 5.5 MeV, m^ = 140.7 MeV, 
and with different flavor-mixing interaction strengthes: i^A^ = 12.36 (thicker curves) and 
i^A^ = 8.0 (thinner curves). Right panel: phase diagrams for m; = m^ = 5.5 MeV (thicker 
curves) and in the chiral limit mi = rUs = (thinner curves), respectively, with KA^ fixed 
at 12.36. 

temperature driven chiral symmetry restoration transition is first order for three massless 
quarks. In the framework of the PNJL model, this could be understood as a result of that the 
thermodynamical potential involves a term cubic in the chiral condensate, because of the 't 
Hooft interaction. To study the the dependence of the phase evolution on the fiavor-mixing 
interaction more thoroughly, we perform a series calculations by reducing the strength of 
the fiavor-mixing interactions gradually and find that when the fiavor-mixing interaction 
strength KA^ decreases to a value of about 5, the order of the temperature driven chiral 
transition is changed from first to second. 

5 Summary and Conclusions 

In summary, we have extended the Polyakov-loop improved NJL model to the 2+1 fiavor case 
with inclusion of strange quark. This 2+1 fiavor PNJL is a synthesis of the conventional 2+1 
fiavor NJL model, which includes the fiavor-mixing 't Hooft interaction, and the Polyakov- 
loop dynamics governed by a Polyakov-loop effective potential. Within the framework of 
such a model, we have studied the chiral and Polyakov-loop dynamics and their mutual 
infiuences to understand the nature of the QCD phase transitions in the three-fiavor system. 
More concretely, we investigate the chiral and deconfinement crossovers with finite tem- 
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perature at zero quark chemical potential with physical current mass of strange quark. Three 
kinds of pseudo-critical temperature corresponding to three different quantities: Polyakov- 
loop, chiral susceptibility of light quarks and strange quark number susceptibility, are de- 
termined in the PNJL model. By employing two Polyakov-loop effective potentials existing 
in literatures, viz. the polynomial effective potential and the improved one, we all find 
that different observables lead to different values of transition temperature due to the non- 
singularity of the crossover, which probably indicate that this phenomenon is independent 
of the choice of the Polyakov-loop effective potential in the PNJL model. However, other 
effective potentials need to be developed to verify this result. The hierarchy in the pseudo- 
critical temperatures found in our model is consistent with the recent Lattice QCD results 
in Ref. [15]. However, Tc{xi), ^c(Xs) ^ind Tc{P) found in another Lattice simulation [25] 
almost coincide and there is not the hierarchy. But we should be careful to use the results 
in Ref. [25], because the lattice spacings used in Ref. [25] are not in the scaling regime and 
the results obtained with their lattice spacings can not give a consistent continuum limit 
for Tc [15]. Making use of the two different Polyakov-loop effective potentials, we find that 
there is always an inflection point in the curve of strange quark number susceptibility vs 
temperature, accompanying the appearance of the deconfinement phase, independent of the 
strength of flavor-mixing interaction, which is also consistent with the results of Lattice QCD 
simulations [15]. Effects of the current mass of strange quark {rris) are studied, too. We find 
that the chiral crossover for light quarks moves down to lower temperature and become more 
abruptly with the decrease of m^, while the deconfinement crossover is almost not influenced 
by the variation of m^. 

Furthermore, predictions for nonzero quark (baryon) chemical potential and finite tem- 
perature are made in this work. We investigate the temperature and chemical potential 
dependence of the Polyakov-loop $, A<1> and the condensates for light quarks and strange 
quarks as well as their mutual interactions. We find that in the deconfinement phase the 
value of $ approaches to one and A$ is much smaller than that of $. We also give the 
phase diagram of the strongly interacting matter in terms of the chemical potential fi and 
temperature T. It shows that the critical endpoint (CEP) moves down to the yU-axis and 
finally disappears with the decrease of the strength of the 't Hooft flavor-mixing interaction. 
On the contrary, the CEP moves toward to the T-axis as the current mass of strange quarks 
is reduced and disappears when the chiral limit is approached. 

It is well known that the Ua{^) symmetry of QCD is explicitly broken by the axial 
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anomaly coming from the gluon dynamics at quantum level. It is very interesting to investi- 
gate the effects of the interplay among the Ua{^) anomaly, the Polyakov-loop dynamics and 
the chiral symmetry breaking at finite temperature and nonzero chemical potential in the 
PNJL model. Following Ref. [35] which studies the properties of pion and sigma mesons in 
the two flavor PNJL model, we can analyze the pseudoscalar mesons and their chiral part- 
ners in three flavor case and their convergence could indicate the restoration of the f/yi(l) 
symmetry with the increase of temperature or chemical potential. These studies are under 
progress and we will report it elsewhere. 
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